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Hydrodynamical SPH simulations of galaxy clusters are used to investigate the metal enrich¬ 
ment of the intracluster medium. The final metallicity abundances of the simulated clusters 
are determined according to the numerical resolution and a number of model parameters. For 
a fiducial set of model prescriptions the results of the simulations indicate iron abundances 
in broad agreement with data. Final X-ray properties are not sensitive to the heating of the 
ICM. This supports a scenario where the ICM evolution of cool clusters is driven by radiative 
cooling. 


1 Introduction 

A large set of observations shows that the ICM of galaxy clusters is rich iri i^tals ( see 10 
and reference cited therein). Abundance gradients have also been measured BoB. For instance, 
Ezawa et al.@ found a decrease in the iron abundance of AWM7 from ~ 0.5 solar in the center 
to ~ 0.2 solar at a distance of ~ 500Apc. These measurements provide strong support for the 
supernova (SN) scenario as heating source for the ICM. 

The relative abundance of the heavy elements can then be used to constrain the enrichment 
mechanism of the ICM and the energy input from SNe. Analysis of the spatial distribution 
of metallicity gradients is also important for discriminating among the proposed enrichment 
scenarios §. Analytical estimates of the energy input can be substantially different according to 
the assumed spatial distribution of the ICM metals and the efficiency of transfering the kinetic 
energy released in a SN explosion to the ambient gas0. Hydrodynamical simulations of galaxy 
clusters have the advantage over analytical methods that they take into account self-consistently 
the dynamical evolution of the gas. The implementation of a metal enrichment model for the 
ICM in hydrodynamical simulations is therefore important in order to investigate the ICM metal 
evolution. Chemical e volut ion in hydrodynamical simulations has already been considered in a 
variety of contexts 






In this paper I present results from hydrodynamical SPH simulations of galaxy clusters, that 
have been used to investigate the dependence of the final iron abundance on a number of model 
parameters that control the ICM metal enrichment. This is in done in order to obtain for the 
simulated clusters a final ICM distribution which can simultaneously fit a set of observational 
constraints, such as the observed iron abundances and the luminosity-temperature relation. 

2 Simulations and results 

Hydrodynamical TREESPH simulations have been performed in physical coordinates for four 
test clusters. A detailed description of the procedure can be found in 111. The cosmological 
model is a flat CDM model, with vacuum energy density Ha = 0.7, matter density parameter 
Hm = 0.3 and Hubble constant h = 0.7 = Ho/100Kmsec~^Mpc~^. = 0.015h~^ is the value 
of the baryonic density. The four test clusters have decreasing virial mass and labels LOO, L39, 
L40 and L119. The virial temperatures range from ~ QKeV (LOO) down to ~ IKeV (L119). 
The simulations have a number of gas particles Ng ~ 22, 600. For L119 high resolution (H) runs 
have Ng ~ 70,000. The cooling function of the gas particles depends also on the gas metallicity, 
and cold gas particles are subject to star formation. Once a star particle is created it will release 
energy into the surrounding gas through SN explosions of type H and la. The feedback energy 
{IQ^^erg) is returned to the nearest neighbor gas particles of the star particle, according to its 
lifetime and IMF0. SN explosions also inject enriched material into the ICM, thus increasing 
its metallicity with time. The adopted SNIa yields are those of Iwamoto et al.0 (model W7). 
The yields of type H SN are those of Woosley & Weaver^ (model B). For each star particle 
the ejected masses of metals are calculated at each time step and distributed among the gas 
neighbors (^ 32) of the star particle. For the smoothing spline a common choice is the SPH 
kernel EjO. The profiles shown in the Figure are for a uniform deposition of metals. 

The projected emission weighted temperature profiles are plotted as a function of the rescaled 
radius in panel (a). Data points are the n^an error-weighted temperature profiles of 11 cooling 
flow clusters from De Grand! &: Molendi Q. Each smoothed profile has been rescaled to match 
the last data point. There is a remarkable agreement of the simulated profiles with data, the only 
important exception being for the two innermost bins. The increase of the cluster temperature 
at the centre is a consequence of the ei^opy conservation during the galaxy formation and the 
subsequent removal of low-entropy gas Eil. A possible explanation for this discrepancy with the 
data lies in the fact that the temperatures being measured are spectral temperatures El. 

The projected metallicity profiles are displayed as a function of distance in panel (b). The 
data points are the mean profile from the nine cooling flow clusters of De Grand! &: Molendi 0. 
The iron profiles of LOO and L39 are the ones which are in better agreement with data. The 
overall shape of L40 is similar to that of the other two runs, but has a lower amplitude. The 
systematic difference between the profiles of runs L39 and L40 is mostly due to the different 
dynamical histories of the two clusters. Therefore, there are uncertainties in the final profiles 
which can be as high as ~ 50% and are related to to the cluster dynamical evolution. At outer 
radii all of the profiles show a radial decay steeper than that of data points, which in fact can 
be considered to have an almost constant profile. It appears very difficult to modify the model 
parameters of the simulations in such a way that the simulated profiles match the data points 
at the outer radii, without also increasing central abundances. 

The iron profile of L119H is inconsistent with the data points, the iron abundances being 
smaller at all radii. This is the coldest of the four test clusters {Tx ~ IKeV). The averaged ob¬ 
served profile is that of 9 cooling flow clusters, with minimum temperatures AKeV. Without 
measured profiles for cool clusters it is therefore difficult to put observational constraints on the 
different model parameters from the simulated profiles. 

Global values Ape = ^Fe/^H of the iron abundances for the simulated clusters (open 





squares) are compared in panel (c) with the estimated values (filled squares) for the nearby 
cluster sample of Matsumoto et al. Ej. There is good agreement with data, but L119H has 
Ape — 0.2, which is about a factor ~ two smaller than that expected from extrapolating the 
sample average below the minimum cluster sample temperature (~ 2KeV ). For cool clusters 
there is the clear tendency in the simulations to produce a smaller amount of iron than that 
inferred from observations. 

Finally, the values of the bolometric X-ray luminosity are shown in panel (d) as a function 
of the cluster temperature. Mass-weighted temperatures have been used as unbiased estimators 
of the sp^tral temperatures 111 . Data points (open squares) are those of Fig. 11 of Tozzi & 
Norman Ea. For the sake of clarity, not all of the points of their Figure are plotted in the panel. 
For a consistent comparison with data, a central region of size b0h~^Kpc has been excised El, 
in order to remove the contribution to Lx of the central cooling flow. The continuous line is 
the best-fit El Lx = 3.11 • 10^^h~^{T/6KeV)‘^'^^. The dashed lines are the 68% confidence 
intervals. The Lx of the simulations are in excellent agreement with the data over the entire 
range of temperatures. An additional run (LllDHb) has been performed for the cluster with the 
lowest temperature, the only difference with respect to L119H being a SN feedback energy of 
IqSO £qj, t)oth SNII and la. As can be seen, the final Lx of the run L119Hb is very similar 
to that of L119H. This demonstrates that final X-ray luminosities of the simulations are not 
sensitive to the amount of SN feedback energy that has heated the ICM. In fact, a simulation 
run with zero SN energy being returned to the ICM yields very amilar results. These findings 
are particularly relevant in connection with the recent proposal H that the X-ray properties of 
the ICM are driven by the efficiency of galaxy formation, rather than by heating due to non- 
gravitational processes. To summarize, the metal enrichment of the ICM can be modeled in 





















hydrodynamical SPH simulations with a number of model parameters. The results show hnal 
iron abundance gradients in broad agreement with the data. Moreover, the final X-ray properties 
of the simulated clusters are not sensitive to the amount of SN energy that has heated the ICM. 
These Endings support a scenario where the ICM evolution of cool clusters is dominated by 
cooling and galaxy formation. 
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